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ABSTRACT 


An analysis and test rig measurements of the per- 
formance of a transonic axial turbine are reported. The 
purpose was to confirm the accuracy of measurements made 
in a test rig which was designed to separate the losses 
occurring in the stator from the losses occurring in the 
rotor blade rows. The analysis was programmed for the 
Hewlett-Packard 21-MX computer. Reasonable agreement 
between predicted and measured characteristics was obtained 
using experimentally determined losses in the computer 
program. Lack of agreement was noted using theoretical 
values. It was concluded that the rotor was not choked 
at the conditions in the tests, and that the test rig 
measurements were valid. A successful technique for 


smoothing the data obtained from the rig is also reported, 





I. INTRODUCTION 


The transonic turbine test rig installation at the 
Naval Postgraduate School Turbopropulsion Laboratory was 
designed to study the effects on turbine performance of 
varying axial and tip clearances, to study the effects of 
blading design on turbine performance, and to allow the 
Separate determination of stator and rotor losses in an 
operating machine. 

Until the work of Solms (Ref. 1) the separation of 
rotor and stator losses thatch test rig measurements had 
not been attained satisfactorily. Through improvements 
in hardware and instrumentation and improvements in the data 
reduction process the stator and rotor losses were deter- 
mined separately and were reported in Ref. 1. 

Anomalies remained, however. Specifically, the turbine 
coifiguration designated as Turbine C in Ref. 1 gave 
different results when compared in terms of "referred" 
quantities depending on whether the discharge was to 
atmospheric pressure or to a region of reduced pressure. 
(Turbine C had converging-diverging stator passages in an 
aXial entry, single impulse stage. The turbine was designed 
to operate in the transonic range.) In addition, Raf. 1 
reported considerable scatter in the loss coefficients. 

In the work of Robbins (Ref. 2) it was shown that 
discharge pressure affected the measurements of flow rate 


into the stage and also affected the labyrinth leak rate. 





Accordingly, Robbins determined accurately the flow rate 
into the stage and the leak rate through the labyrinth seal 
(See Fig. 1) for all operating conditions. . The results | 
then obtained for turbine C were reported fully in Ref. 2. 
The continued presence of scatter in the measured loss 
coefficients was reported and a smoothing technique to 
eliminate the scatter was suggested. 

Before the Turbine Test Rig could be used to measure 
the effects of varying parameters: 

1. The scatter in the loss coefficients had to be 

eliminated. 

2. The overall accuracy of the performance results 
evaluated from the rig measurements had to be 
verified in some way. 

The resolution of these problems was the goal of the 

present work and is the subject of this report. First, 
a satisfactory method was found for smoothing the loss 
coefficients. The method is described in Section III. 

The approach taken to verify the performance of the 
rig was to first devise an analysis which predicted the 
performance of the turbine in terms of unknown loss coeffic- 
ients, and then to show that the measured loss coefficients 
were consistent with the predicted behaviour. 

A description of the Turbine Test Rig is given in 
section II. The analysis of the behaviour of the turbine, 
involving a comparison of an analytical prediction with 


the results of a short test program, is described in 
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Section IV. Details of the analysis and the computer 


program are given in Appendix A. 
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Il. TURBINE TEST RIG INSTALLATION 


A. DESCRIPTION 

The test installation consists of three major compon- 
ents; an Allis-Chalmers twelve stage axial flow compressor, 
an exhauster assembly, and the turbine test rig (TTR) 
itself. 

The compressor is the source of driving air for the 
TTR and for the exhauster assembly. Fig. 2 shows the 
piping arrangement. Turbine air passes through the first 
settling tank into an eight-inch pipe containing a flow 
nozzle, into the second settling tank and into the turbine. 

Fig. 3 shows the plenum, the floating stator assembly, 
the rotor, and the dynamometer (Ref. 1). Pressure ratios 
of 6:1 can be achieved when the system is hooded. The 
hood was needed to achieve high pressure ratios in the 
tests reported here. Fig. 4 shows the turbine blading 
of the stator and rotor. Ref. 3 contains detailed descrip- 
tions of the test rig hardware. 

The floating stator assembly shown in Fig. 3 permits 
measurements of the axial force and the torque on the 
assembly. Axial and rotational movements are constrained 
by calibrated force transducers that are heat insensitive. 
These measurements, together with wall static pressure 
measurements, allow the determination of the average 
axial and tangential velocity components at the stator exit. 


In this report one configuration designated Turbine C 
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was tested, the geometry of which is shown in-Fig. 5. 

Table I describes the geometry quantitatively. The stator 
blade prof#le is shown in Fig. 6. The blades of the stator 
generate a converging-diverging nozzle shape. Pressure 
measurements were taken at the locations shown in Fig. 6. 
The pressures necessary to the analysis of the stator axial 


force were taken at the locations Shown in Fig. 7. 


B. TEST MEASUREMENTS AND ACCURACY 

1. Mass Flow Rates 

Appendix A of Ref. 2 gives a detailed description 

of the method used to determine both the turbine flow rate 
and the labyrinth seal leak rate. 

2. Forces, Torques, Temperatures, and Pressure 

Ref. 3 and Ref. 5 give calibration procedures 

for the TTR. Identical procedures were employed here. 
Table II of Ref. 1 gives the expected accuracies of the 


measurements. 


C. TESTING AND DATA REDUCTION 

The TTR data collection system is described in Ref. 4. 
Appendix D of Ref. 1 gives a detailed explanation of the 
turbine test procedures. Those procedures were followed 
here with the exception that a constant RPM was held and 
the pressure ratio varied over the desired range. The 
data reduction method developed in Ref. i and Ref. 2 was 


revised as described in Appendix B. 
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Jit. DATA SMOOTHING TECHNIQUE 


The sensitivity of the loss coefficients to variations 
in measured quantities was shown in Ref. 2. It was also 
stated that some variation is unavoidable since measure-~ 
ments are taken over a period of more than a minute. It 
was also pointed out in Ref. 2 that the parameter most 
important to the calculation of the loss coefficients was 
P, (the average pressure at the stator exit). P, is an 
average pressure that can not be measured directly. It 
is derived as aegeniped in Ref. 1 from many other measure- 
ments. Significant scatter was observed in the variation 
of P, as speed was varied at fixed pressure ratio. However, 
it was found that the hub and tip pressures (Py and P.) 
measured just downstream of the stator varied smoothly over 
the same range. Since these two pressure were measured 
directly, and since they varied smoothly, it was assumed 
that the pressure behind the stator must vary smoothly also. 
Consequently, it was determined that a polynomial curve fit 
could be used to describe the variation of P, aS speed was 
varied. The variation of P, was represented as a function 
of P,, P,, and the isentropic head coefficient (K, 5) 
in the forms 


GC - Rape = AL tALKis +A, ee (1) 


Where A; is the polynomial coefficient. 


Then, oe (where P,, is total pressure upstream of 


the stator) was computed using the expression 
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or =i PtTth Peyeeie 
Peo Pto P50 

Fig. 8 is an illustration of this procedure, which 
was added to the "bulk process" data reduction program. 
The data points in Fig. 8 were taken from Runs 6 and 7 
in Ref. 2. The two Runs were for the same poncietens but 
were made at different times. It can be seen that the 
trends are the same for both runs but the scatter is 
considerable. The scatter is the cause of the scatter in the 
calculated loss coefficients. The lines on Fig. 8 are the 
polynomial approximations according to Eq. (1) for the two 
runs. It can be seen that the polynomial approximation 
averages the data and maintains the original trend. The 
chosen smoothing function was a ee degree polynomial. 

Fig. 9 shows the stator loss coefficient prior to 
smoothing for the points in Run #7. Fig. 10 shows the 
Same data after smoothing. 

Similarly, Fig. 11 shows the rotor loss coefficients 
for Run #7 prior to smoothing and Fig. 12 shows the same 
data after smoothing. 

It can be seen from these figures that the scatter 
has been removed. This is particularly true for the rotor 
loss coefficients. As pointed out above, it is believed 
that the observed scatter was due to the sensitivity of 
the loss coefficients to small changes in measured quan- 
tities during the data collection process. The smoothing 


technique removes the random variations recorded during 
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data collection and results in a much more realistic 
representation of the variation in the losses. 
The smoothing technique was incorporated into the 


data reduction program as described in Appendix B. 
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IV. ANALYSIS OF TURBINE PERFORMANCE 


A. APPROACH 

In order to determine if the performance evaluated 
from the rig measurements was accurate an analysis to 
predict the behaviour of the test turbine was carried out 
and programmed in BASIC language. 

A performance test was then conducted in a particular 
way in order to provide a pamaeesaer of the measured with 


the predicted behaviour. 


Bs ANALYTICAL PREDICTIONS 

The flow through the turbine was analysed using a 
pseudo-1 Dimensional compressible approach. The analysis 
is described in detail in Appendix A, together with the 
computer program which was used to obtain predictions of 
the turbine performance. One of the inputs which the 
program requires is the rotor passage loss coefficient at 
zero incidence. Using the method given by Vavra in Ref. 6 
it was determined that the rotor loss coefficient should 
have the value .2514%. However, results of previous tests 
of the turbine indicated that the rotor loss coefficient 
was rarely as high as .2514 and could be as low as .1. 
Therefore, the performance of the turbine was analysed 
using rotor loss coefficients of .2514 and .1. The pre- 
diction program was run for both 15,000 and 18,000 RPM with 
an assumed rotor loss coefficient of .2514, and for 18,000 


RPM with an assumed rotor loss coefficient of .1 
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The above parameters were chosen to obtain a prediction 
of the turbine performance in a range in which experimental 
data could be obtained. In particular, the analysis could 
be used to predict ne pressure ratios at which choking 
would occur in the rotor as well as in the stator. The 
pressure ratio at which choking occurred in the rotor could 
then be established experimentally in a test conducted at 
fixed speed. An examination of the choking condition was 
considered to be a first test of the performance analysis. 

The results of the analysis for the parmameters given 
above are shown in Fig. 13, Fig. 14, and Fig. 15. What is 
shown in the figures is the map of the range of values 
which unknown parameters in the analysis can have, that 
leads to a solution. 

Figures 16-19 show predicted performance parameters 
for the case of 18,000 RPM and assumed rotor loss of .1. 
These results will be discussed in conjunction with the 


results of the turbine test run. 


C. EXPERIMENTAL TURBINE TEST 

In order to examine the occurrence of rotor choking 
the turbine was run at constant RPM and the pressure ratio 
across the stage was increased in increments by lowering 
the back pressure. The point at which the horsepower ceased 
to increase for an increase in pressure ratio was examined 
to determine the choking point. At the condition where 
the flow reaches a Mach number of unity at the exit of the 


rotor, the power produced by the turbine can not be changed 
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by altering the downstream pressure. 

Tests were conducted in this manner at 15,000 and 
18,000 RPM. The controlled parameters for the tests are 
given in Table II. The reduced data is given in Table III. 
The referred horsepower is shown plotted versus the pressure 
ratio for the 15,000 RPM run in Fig. 20 and for the 18,000 
RPM run in Fig. 21. 


The results are discussed in the next section. 
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D. COMPARISON OF ACTUAL AND PREDICTED PERFORMANCE 

Figures 13, 14, and 15 are maps of possible solutions 
for the flow through the turbine, when the parameters which 
are unknown are allowed to vary. It is noted first that 
the predicted range of solutions extends to pressure _ 
ratios Cyl below the predicted choking line. The 
explanation for this apparent anomaly is that the program 
calculates all possible solutions, and for any point below 
the line a particular combination of losses and blockage 
factors existed which would allow a solution at that point 
without choking. It should be pointed out that any point on 
the plot can be brought to the choking line by reducing 
the stator loss by a very small amount. This is the 
procedure that was followed to get the range of values 
at choking shown in Figures 16-19. 

Fig. 13 is the map produced by the program for RPM 
equal 15,000 and for an assumed rotor passage loss coeffic- 
jent equal to .2514. As can be seen from the figure the 
pressure ratio, (P,/P,,)> at choking was about .28, cor- 
responding to a stage pressure ratio (Bee/2e)) of 355517 
Fig. 14 is the map for RPM=18,000. Note that the predicted 
Stage pressure ratio for choking is again about 3.57. 

Fig. 20 and Fig. 21 show the variation of the referred 
horsepower vs. the pressure ratio which were measured in 
turbine tests at 15,000 and 18,000 RPM respectively. It 
can be seen that the referred horsepower did not become 


independent of the back pressure at the pressure ratios 
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achieved in these tests. It was concluded that the stage 
did not choke at the predicted pressure ratio of 3.57. 

It was noted however that at the highest pressure ratios 
obtained at 18,000 RPM, the slope of the horsepower curve 
was becoming smaller. 

The velocity diagrams in Figures 23-31 are also consis- 
tent with the argument that the rotor did not choke during 
the turbine tests. The velocity from the rotor, Vis 
behaves smoothly and ina predictable manner throughout 
the pressure range at both test speeds. This might not 
be expected if, following choking, shock waves appeared 
downstream of the rotor exit plane. 

Note the values of rotor loss coefficient given for 
the test results in Table III. With the exception of 
point #10, which was considered to be in error, they were 
all smaller than the value (.2514) assumed in the first 
performance calculations. As explained in Appendix A, 
the value of the rotor passage loss coefficient used in 
the program is the smallest value that can be calculated 
for the overall rotor loss coefficient. Therefore, if 
the experimental results were accurate, the computer program 
could not predict the correct performance since it could 
never calculate a rotor loss less than the input value of 
the passage loss, which was .2514. It was as a consequence 
of this observation that the performance was re-calculated 
using a rotor passage loss coefficient equal to .1. 


Figures 16-19 show the choking behaviour of the 
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turbine predicted using the computer program for an assumed 
rotor passage loss of .1. In these figures, the data for 
test point #6 is also shown. In this case it can be seen 
that the predicted pressure ratio era) at choking is 
about .245, corresponding to a stage pressure ratio (EEe EE) 
of about 4.08. The latter is slightly higher than the 
maximum pressure ratio that was attained in the turbine 
tests. (At 18,000 RPM the highest attainable pressure 

ratio was 3.97, due to the characteristics of the dynamom- 
eter.) 

Fig. 17 is the range of efficiencies predicted for the 
choked condition. It is noted that the efficiency measured 
at test point #6 was reasonably close to the predicted 
maximum efficiency. 

Fig. 18 shows the range of stator losses for which 
solutions existed, in comparison with the value of stator 
loss measured at test point #6. It can be seen that the 
measured value intersects the predicted range of possible 
solutions. Fig. 19 shows the predicted range of the rotor 
coefficient. Again, it can be seen that the measured value 
at test point #6 overlaps the predicted range of possible 
solutions. 

Fig. 16 shows the predicted range of horsepower 
compared with the horsepower measured at test point #6. 

The measured horsepower was Slightly greater than the max- 


imum value which was predicted. 
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ye DISCUSSION 

The results shown in Figures 16-19 illustrate the 
uncertainty in the prediction of the performance of the 
turbine when only the rotor passage loss and blockage 
factor are known. It is recalled that the analysis 
Satisfies only continuity through the stage, and values of 
an additional loss coefficient and an additional blockage 
factor must be established to obtain a unique solution. 

It is the aerodynamic shaping of the surfaces which deter- 
mines these factors. Figures 16-19 show the possible 
range of performance that results simply from the areas 

of the passages and the blade angles. 

The irregular shapes of the bounds on the possible 
solutions in Figures 16-19 are interesting. There is no 
obvious explanation for the reduced range of solutions near 
Kut = ,81. 

The results obtained in the present work suggest 
that the performance of the turbine as measured in the 
turbine test rig is reliable. It had been thought prev- 
iously that the rotor loss measurements were too low. 

Here, an analysis was carried out to predict the performance 
of the turbine, and the predicted results have shown very 
good agreement with the performance measurements in the test 
me 

Whether or not the rotor was choking had also been a 
question in the past. In the present work both the analysis 


and experimental tests have shown that the rotor was not 
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choked in any test to date. The computer program with an 
assumed rotor loss of .1 has shown that the test rig should 
not be choked at any pressure ratio below about 4.08. Test 
data has shown that the rotor was not choked at a pressure 
ratio just slightly below 4.0. 

It is also of interest to note that for a rotor 
passage loss coefficient of .2514 the predicted pressure 
ratio was shown to be independent of RPM. At both 15,000 
and 18,000 RPM the choking pressure ratio (a= / = was 
calculated to be about 3.57. More results are needed to 
confirm that the choking pressure ratio is indeed indepen- 


dent of speed. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The rotor of Turbine C was not choked at any pressure 
ratio tested so far. This conclusion is based on the 
prediction of the computer program and also on the 


turbine test results. 


The loss coefficients measured in the TTR and smoothed 
as described in this report can be accepted as being 
truly representative of the losses in the separate 


blade rows of the turbine. 


In particular, the results from the computer program 
Suggest that the magnitude of the measured rotor loss 
coefficients is probably correct. The method used 

by Vavra in Ref. 6 predicts much larger values of loss 
coefficients for the present rotor geometry than those 
which were measured. When the measured value of the 
rotor loss coefficient was entered into the computer 
program, there was good agreement between all the 
calculated and the measured turbine performance 
parameters. When the higher (calculated) value of the 
rotor loss was entered, there was a pronounced dis- 
agreement between the calculated and the measured 


performance. 


The computer program shouid be used (and progressively 
developed) in conjunction with all tests carried out 


in the turbine test rig. Most importantly, a test 
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should be conducted to determine the rotor choking 
condition experimentally. This will require the purchase 
of a water-brake dynamometer to extend the power-speed 


range of the test rig. 


Since the accuracy of the test rig measurements is 
no longer in doubt, experiments to determine the effect 
of parameter changes (e.g. axial and tip clearances) 


on turbine performance can go ahead. 
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TABLE I 


TURBINE GEOMETRY 





a. | | a: ae 
DESCRIPTION | aA* ; aoa Ry 
(in) | “BLapEesh  _ (IN) 


CONVERGING 
STATOR 1 | DIVERGING 2.9058 34 4.184 10.5775 
NOZZLES 
h 0.732 
740 ns 
=0.8475 
m 






















2? 





eA untae 
PARAMETERS FOR TURBINE TEST 


por @ Fes 


1 15,000 2.01 
2 1 5), O00 2.42 
3 Ton Oce Zao 
4 ow 3.49 
5 15,000 3.68 
6 13, COO 3.94 
7 18,000 3.47 
8 18,000 2.98 
9 18,000 ee 
10 18,000 Zoe 
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FIGURE 10 STATOR LOSS VS. ISENTROPIC HEAD 
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FIGURE 12 ROTOR LOSS VS. ISENTROPIC HEAD 
COEFFICIENT, AFTER SMOOTHING 
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FIGURE 17 PREDICTED EFFICIENCY Vs. Kod 


| [BA 
Aba Baa 98° ha’ co a Aaa BL°A 9L°A hL° aL ALR 












So°B 
; + T°=SSOT YOLON CaNNSSV 
000 ‘ gT=Wd¥ 
tYOd GELOIdgaYd 
; 1 Zt °DId 
: 69°2 
+ + 
+ 
4 + 
EL‘ 
+ 
+ 
7 LL’@ 
t+ + ¢ + $F $F F SF + t+ 4 
1a°2 


PREDICTED EFF IECIENCY 
50 





FIGURE 18 PREDICTED STATOR LOSS VS. Kat 
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FIGURE 19 PREDICTED ROTOR LOSS Vs. Ryt 


= 
+ 
= 
po 
— 
ae 
——_ 
_ 
=> 
+ “4 
+ + 
+ + 
+ + oc 
© 
“ © 
a 
— ON 
ei 
SS he) 
a ee 
- 
o fy 
a % 
; i a z 
BR] BJ BJ BJ a | BJ 


PREDICTED ROTOR LOSS COEFFICIENT 


by 


RPM=18 , 000 


=.1 


ASSUMED ROTOR LOSS 





0.72 8.74 8.75 2.78 2.88 8.82 2.84 8.86 2.88 4.50 
KB | 


4.70 





FIGURE 20 REFERRED HORSEPOWER VS. PRESSURE RATIO 
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FIGURE 21 REFERRED HORSEPOWER VS. PRESSURE RATIO 


3.408 3.20 3.448 3.68 3.88 4.02 


PTO/P2 


2.80 





2.50 





2.48 





2.24 


21 
REFERRED HORSEPOWER 


PRESSURE RATIO 


FOR RPM=18 , 000 






FIG. 





VS. 
2.00 


. 82 


 ) 


SJ 5s 


BJ 
sx ™ 


MEASURED REFERRED HORSEPOWER 


54.02 
2 
20. He 


54 





FIGURE 22 VELOCITY DIAGRAM FOR TURBINE TEST 
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FIGURE 24 VELOCITY DIAGRAM FOR TURBINE TEST 
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FIGURE 25 VELOCITY DIAGRAM FOR TURBINE TEST 
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FIGURE 26 VELOCITY DIAGRAM FOR TURBINE TEST 
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FIGURE 27 VELOCITY DIAGRAM FOR TURBINE TEST 
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FIGURE 28 VELOCITY DIAGRAM FOR TURBINE TEST 
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FIGURE 29 VELOCITY DIAGRAM FOR TURBINE TEST 
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FIGURE 30 VELOCITY DIAGRAM FOR TURBINE TEST 
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BUG URES 31. 
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APPENDIX A 
TURBINE PERFORMANCE ANALYSIS 


A-1 INTRODUCTION 

This appendix gives a detailed description of the 
theory and analytical technique used to predict the 
performance of a turbine in the Turbine Test Rig. Variables 
used in this appendix are defined as they are introduced. 
Also given are the details of the computer program used 
to implement the analytical method. The variables used 


in the computer program are listed and defined in Table A-I. 


A-2 ASSUMPTIONS 

The following assumptions were made: 

1. The stage was choked in the stator. 

2. In solving the continuity equation between the 
stator and rotor the solution corresponding to the lower 
Static pressure was taken. 

Assumption (1) was made on the basis of past exper- 
ience with the test rig. Ref. 1 states that there is 
Supersonic flow in the stator and an examination of the 
pressure distribution in the stator passages confirms this 
statement. Fig. A-1 is a plot of the pressure distri- 
bution through the stator by pressure tap number (See Fig. 6). 
It can be seen in this plot that there was a sharp pressure 
rise between tap #3 and tap #4, which indicates the 
presence of a shock in the divergent section of the passage. 


Also, the level of static pressure at the throat taps 
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in comparison to the supply pressure, was observed to be 
independent of downstream pressure for all operating con- 
ditions. 

Assumption (2) limits the range of possible solutions 
as shown in Fig. A-2. This figure is a plot of flow 
coefficient (%) vs. pressure ratio CA for various 
loss coefficients (z). (The flow coefficient, a non-dimen- 
Sional mass flux, is introduced in Section A-2.) There 
are two possible pressures for a given flow coefficient 
and given loss coefficient. Solutions were limited to 
the left side of the line connecting the points of maximum 
flow function. The procedure was chosen because it resulted 
in the “supersonic root" for the flow from the stator. 

Also, conditions corresponding to choking were sought 

for the rotor, and the locus of locally sonic conditions 
is aS shown in Fig. A-3. It should be noted that because 
of the definitions of the flow coefficient and loss coef- 
ficient, the "choking" condition at which the downstream 
pressure has no effect upstream of the plane in question, 
does not correspond to the maximum mass flux from given 
stagnation conditions. 

The analysis is pseudo-1i-dimensional. Blockage factors 
are introduced with the physical cross sectional area to 
account for non-uniform flow conditions. Loss coefficients 
are defined on the basis of kinetic energies. The analysis, 
which requires only mass flow continuity through the stage, 


is divided into three stages; the stator, the interblade 
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region and the rotor. 

While the analysis itself is general, solutions can 
only be obtained by assigning values to chosen parameters. 
In the solutions reported here, the flow angle at the stator 
exit (4;), the blockage factors for the stator throat | 
(ke. ) and rotor exit plane (Kip) s and the rotor passage 
loss coefficient (z,) were given values. The value of Xy 
was fixed atay =75° based on the results of the turbine 
tests and from the blading geometry. The blockage factors 
were set at Kk o= +965 and Ky p= 085. The rotor passage loss 
coefficient (zp) was an input variable which was set at 
22514 (theoretical value) or .1 (lowest value obtained in 
test results). The "rotor loss coefficient" output from 
the program, Za oe includes all losses from the stator exit 
plane to the rotor exit plane. 24.9 corresponds to the 
"rotor loss coefficient" evaluated from test rig measure- 


ments: consequently, 24.2> Zp 


67 





A-3 BASIC RELATIONS 


The sketch 


aA 





shows a general 


adiabatic process 
with entropy increase 
(losses) on a T-S ie 


; : The 
diagram on which Zs 


Xis 
the temperature = | 


scale has been divided S 


() 


wo 
San sewage) 


Ny 


2 
tS 


ca 


by the constant stagnation temperature of the process. A 
perfect gas is assumed. (0) represents the initial stag- 
nation state and ( ) represents the final state. 

The non-dimensional velocity, X, is defined as 
X “yh where V,. =v2C,T,is the "total" or "limiting" 
velocity, and the loss coefficient, z, is defined as shown. 
Note that the stagnation velocity is constant throughout 
the process. 

The "flow coefficient", %, is defined here as 


Ue a 
Oe TAWA, = 


where Ww is the flow rate, (is the density at the initial 
stagnation temperature and pressure, A is the flow area and 
Ky is the blockage factor, The flow coefficient defined 
in this way is a non-dimensional flow rate per unit area, or 
mass flux, referred to initial stagnation conditions. 

The flow rate is given by 


wa PAR V 
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so that Eq. A(1i) can be written 


7 “4 Veo con 


which, in terms of the non-dimensional velocity defined 


above, becomes 


al i | 
ee) (ae - 


The loss coefficient is defined, as shown in the above 
Sketch, as 
pee ic ea 
n= = /— : A€3) 
Wey -~Tis Xi 
Rearranging Eq. A(3), 


2 2 
(-2)* = Ee 
a Tee 
and using the isentropic relationship between pressure and 


temperatures 


= (0 2)[1- &)* | ACH) 


Using Eq. A(4), Eq. A(2) becomes a single equation 
for @ as a function of pressure ratio DAD In analysing 
the turbine, generally the flow rate is known so that the 
flow coefficient can be calculated using Eq. A(1). On 
specifying a value for the loss coefficient,’(z), the pres- 
sure ratio can be calculated from Eq. A(2) using Eq. A€4). 
An iterative technique is used (Newton's Method) starting 
with an initial estimate of the pressure ratio. Fig. A(2) 


shows @ as a function of Aan for various values of z, 
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in graphical form. The solution to the left of the maxima 
is obtained by beginning the iteration with a small value 
of P,/P,.- 

The analysis of the turbine begins by specifying that 
the flow rate is set by the choking of the stator nozzles. 
The pressures at successive stations are then calculated 


in turn from the flow coefficient as described above. 


A-4 STATOR EXIT CONDITIONS 

The flow coefficient at the throat of the stator is the 
value of the flow coefficient corresponding to locally 
sonic conditions at the minimum area. In the present work, 
the maximum flow coefficient at a loss coefficient equal to 
.05 (defined as S*) was taken at the stator throat, 


throughout the calculations. 


Since 
fe = 2 A(5) 
Ses Veo As hy s 


where A. and Kins are the area and blockage factor at the 
throat station, the flow rate being constant through the 


machine requires that 


W = ee h,. O* = constant 


at stations ahead of the rotor, where T 7 1S constant; this 


¢ 
implies that 
ee (ie hips * = constant 
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The value of £* can be obtained analytically when the 
loss coefficient is known. The maximum occurs where 
ae O, which leads to the solution (denoting values at the 
aOR throat by an asterisk), 


(x*)* = B -yB* -.¢ A(6) 
vnere p= 14 GED AG =e -z) (2) 
ne 


J~/ 
and P*/P,. = io | A(7) 


where D = (ar (12S 


Then £* is given by Eq. A(2). 
With £* known, continuity of flow rate requires that 


at the stator exit, 


B, = O* cia Sea A(8) 


lah Ry, Cog a, 


where subcript 1 represents the stator exit plane (station 1) 
and &,is the flow angle at thestator exit. 

Knowing Bs ee can be found by iteration as 
described in section A-3. Solutions on the left side of 
Fig. A-2 were selected by beginning the iteration with 
BAP +o close to zero. 


After obtaining thevalue for Eaten then X, can be 


1 


ae 





calculated using 





x, A(9) 
Then, P,/P +4 can be found from 
oven 
Eyes (teeny A(10) 
and 
eee E/E ee) (eae) A(11) 


A-5 INTERBLADE SPACE 
The interblade 
Space, from Station 


ieee co station (i), 





is shown in the sketch. STATOR 


Conditions at (i) 
are calculated ina manner 


similar to those at Station (1). 


First, by continuity, 


L 


"inp 


6. = ge Ache LA A(12) 
: Ai k,.Cosdd Pe, 


where OA is the flow angle at the rotor entrance. 


Ve 





The value of &{ is obtained from the condition that ang- 


ular momentum is conserved in the interblade space. Then 


R 
X. = eh sink, A(13) 


Ul 
a 


where R,=mean radius at Station 1, and R,=mean rag .useau 
Swacilon i. 


Then 
ee 
eee oul A(14) 
and, also 


= _(IN_.p. 
i - (EH RK, A(15) 


where N is the RPM, and the tangential velocity (Via) 





has been non-dimensionalized as 
X , =u A(16) 
Knowing B., a value for Eee can be found if a 
value of the interblade loss coefficient (z, ) is assumed. 


Then 


Pa / Py = Coa (Eee) A(17) 
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and 





A(18) 


| 


With the absolute flow properties established at 


BVT 
Pe 


the rotor entrance, the relative flow conditions can be 
calculated. First, from the geometry of the velocity 


diagram, the relative flow angle (@;) is given by 





@, = tan? |Acswai-Ui =] A(19) 
x Cos Ac 
where 
he 
GU. =—t A(20) 
: Vig 


is the non-dimensional rotor speed at radius Re. Then the 
non-dimensional relative velocity (Xe) is, from the vel- 


ocity diagram, given by 


_ cos rc 
Wis cos eC: A(21) 


and the temperature Tes by 





A(22) 
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The equivalent temperature (Ref. 1), is given by 


2 
Ton = R 
a= i- x¢ a5 x*, + US (2) - | A(23) 
T i wi ee 
to 1 
and therefore 
¥ 
ceey Pb a ee T Pp al 
to 1 to to “i to 
where 
To; = equivalent temperature into rotor (See Fig. A-7) 
Po, = equivalent pressure into the rotor (See Fig. A-?) 


A-6 ROTOR CONDITIONS 
Using the values calculated thus far it is now possible 
to calculate the conditions in the rotor. Continuity 


requires that 


go, = ge te h A(25) 
TE Cl 


where a subscript R denotes the exit plane in the rotor 
blading. 

It should now be noted that the process in the rotor 
frame from the equivalent conditions at (i) to the exit 
of the rotor is entirely similar to the basic process 
described in Section A-3. However, the equivalent temper- 
ature takes the place of stagnation temperature and 


relative velocity replaces velocity in the given equations. 


- 





The pressure ratio which satisfies the flow coefficient 
given by Eq. A(25), for an assumed rotor passage loss 
coefficient (z,) is Po/P a; This is obtained as described 
in Section A-3. The corresponding non-dimensional velocity 


is now 


A( 26) 





where 


Y,-=— ICAP ) 


Hence 


Ww w2 Pt 


From the velocity diagram, 


x5 = (U, - X,, sin@)* + (X,, cos@,)* A(29) 


where @, = 71° (from the geometry of the rotor), and 





ee tT. ce . 
2 ae R, 
The temperature at the rotor exit is given by 
es 
iE Je y Tae 
ra = {(1-z,) € + Zp ry A(30) 
to BE} to 
a iW 
we a 2 
that mm = mM + X A(31) 
SO a i ite 2 
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Next the loss coefficient that includes all of the 
losses in the stage from Station 1 to Station 2 can be 
calculated. This coefficient includes all of the losses in 
the inter-blade reson as well as those in the rotor 
passage. 1% is called Zs 9 here but, in fact corresponds 
to the measured rotor loss coefficient in the TTR since all 
losses aft of the stator in the TTR are included in the 
"rotor loss coefficient." 


TET 


= a seen 
1-2 i -T, 





| 


Wegener e iegee 
ltée te 
Se te (32) 
= A 
Te (BR sfs\ 7 
Jee ae Pe; 


A-7 STAGE PERFORMANCE PARAMETERS 


il 
+ 


In the notation used for the turbine test rig in 
Ref. 1 and Ref. 2, the following equations determine the 
turbine performance parameters from quantities calculated 


in the above steps: 








T 
+2 
an, ~ Mo (1 ) A(33) 
-/ 
Fe are | 
ee ele ae A (34) 
T 
Nog * a A(35) 
os is 


aa 





Ss 
UH} 


oe (S toV toto kyo) 


= ¥ 8 
H.P. = w (.2h02) ane 





M = H.P, 0) 360. 
(55 aN 
Af Seto. 
Pref 
_ To 
Tet 
we = wie 
5 
H.P.* = HPs 
J fo 
Me = 
a — vee 
Je 
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A(36) 


A(37) 
A(38) 
A(39) | 
A(40) 
A(41) 
A(42) 
A(43) 


A( 44) 





A-8 THE COMPUTER PROGRAM 

A-8.1 Overview 

There are nine degrees of freedom in the com- 

puter program: 

1. A, , the flow angle Outlot the Stator 

Ze Koa the stator exit blockage factor 

aie Kiss the interstage blockage factor 

; KR? the rotor passage blockage factor 


eZ. the stator loss coefficient 


os? 


os 
5 
6% z., the interstage loss coefficient 
ie Zp» the rotor passage loss coefficient 
8. N, turbine RPM 
9. Ee total temperature upstream of the stator 
%, is input on the basis of turbine test results. 
At high pressure ratios the value of A, was measured to 
be US =76°. Therefore, the value of &%, was set and held 


at nS°.. The values of RPM and T are thot conditions 


to 
for which a solution is being sought. 

All the other parameters, the rotor passage loss 
coefficient (zp) and blockage factor (kip) were given 
chosen values. This choice was made because the incidence 
losses were included in the interblade calculation. It 
was thought that the parameters describing the flow inside 
the rotor passages could be held fixed as the speed 
was allowed to vary. 


Consequently, there are four degrees of freedom with 


which the program works. The goal for the program is to 
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find those solutions which result in flow coefficients in 
the rotor that are less than, or equal to, the "choking" 
value. This is done in the program by iteration on the 
remaining four degrees of freedom. Fig. A-5 is a block 
diagram showing schematically the iteration process. The 
order of iteration is (from the most frequent to the least 


frequent), z:, Z 


; 5? Kya? and Kas The loss coefficients for 


both the stator and interblade always begin at zero and they 
are incremented as the program seeks solutions. A graphical 
representation of a typical iteration is shown in Fig. A-4. 
It can be seen that the program can find a solution for 

the condition where the loss coefficient is equal to zero 
(denoted as (1) on Fig. A-4). On the next iteration the 
program will set the loss coefficient equal to .1 and it 

can be seen that a solution still exists at point (2). On 
the next iteration, however, it can be seen that a solution 
does not exist for a loss coefficient equal to .2. In this 
case the program will go back and add a smaller increment 

to the previous loss coefficient (.1). 

Each time an added increment results ina loss coef- 
ficient greater than the maximum allowable the program 
subtracts that increment and adds a smaller one. In this 
way the program eventually finds the maximum loss coef- 
ficient which will yield a solution. At this point, the given 
flow coefficient is at the maximum point for the calculated 
loss coefficient line. At this point "choking", as defined 


here, has occurred. 
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It should be pointed out that this "choking" con- 
dition (the point of maximum flow coefficient for any 
given loss coefficient) is not the point corresponding 
to locally sonic conditions. In Fig. A-3 the locus of 
points of maximum flow coefficient and the locus of points 
of local sonic flow are both shown. The point of locally 
sonic flow always falls to the left of the "choking"™ point 
(as defined in this work), and the flow coefficient for 
sonic flow is always less than the flow coefficient at 
the defined choking point. The reason is that the point 
of maximum mass flux (maximum flow coefficient) is for 
given upstream stagnation conditions (see next section), 
Because of the way the flow coefficient is defined, the 
maximum for a given loss coefficient will not correspond 
to local sonic conditions wxcept where the loss coefficient 
is equal to zero. 

The difference between the "choking" flow coefficient 
corresponding to the maximum and the flow coefficient at 
the local sonic condition is about 10% at a loss coefficient 
of .3 (the highest rotor loss coefficient predicted by 
the program). The difference is about 4% at a loss coef- 
ficient equal to .25 (the highest stator loss predicted), 
and insignificant (less than 1%) at a loss coefficient of 


aieout .15. 
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A-8.2 Description of the Program 
The following discussion refers to Fig. A-6, 
which is a detailed flow chart of the program, to 
Table A-I, which is a list of program variables and their 
definitions, and to the program listing given in Section A-10. 
The first step in the program is to input the variables. 
O&, is input in line 280 (it is also converted to radians 
in that line). k., is input in line 286 and ky 


bi 
in line 284, All other variables are input as requested 


: is input 


by the program (lines 20-110). Note that P, and P, are 


2 to 
input also. They are not involved in the calculations and 
were originally input as the conditions for which a 
particular Solution was sought. 

The next step is to go to a subroutine to calculate 
the value of @* (line 170). £* is defined as the value 
of the flow coefficient at the throat of the stator which 
sets the values of the flow coefficients at all downstream 
stations. §* is generated in the subroutine by the method 
covered in Section A-4 assuming a loss coefficient of .05 
to the throat of the stator (line 160). 

Next the program calculates the conditions in the 
stator exit plane. This is done by first calculating the 
value for ~, (line 320). As explained earlier, the value 
for B, determines the upper limit of the stator loss 
coefficient. 

Using the method described in Section A-3, the program 


calculates the value for P,/P.o° Note that if the stator 
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loss coefficient is large enough, so that a solution does 
not exist, this will be detected in the decision, "ea 
or>1" (lines 350 and 370). If “yes" then the stator loss 
(z,) is too large and the following steps are taken: 

1. The last increment added to the stator loss is 
subtracted (line 420). (Note, the initial stator loss 
increment is .1.) 

2. The existing stator loss increment is multiplied 
by some factor less than 1 (line 430). 

3. The new, Smaller, increment is then added back 
to the stator loss coefficient (line 1570). 

If stator loss increment is sufficiently small 
(line 1575) then the upper limit of the loss coefficient 
has been reached as discussed earlier. At this point no 


further solutions are attainable and Ky is incremented 


a 
(line 2290). Note that the upper limit of stator loss 
will not normally be attained until several iterations 
have been made through the entire program. 

Next the interblade calculations are made. The sol- 
ution for es found for the stator exit plane is used 
to make the calculations for p. (line 770). Again, p. 
determines the upper limit on the interblade loss. In 
the interblade region, however, there is another criterion 
that must be met. The value for 6: generated within the 
program must be equal to @, (8, =69° and is fixed by the 


geometry). This criterion is met by iterating the inter- 


blade loss coefficient as follows: On every pass through 
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this section the interblade loss is initially set to zero 
and the interblade loss SACROTeAt is inttialivese: come 
(lines 720, 740). The value for P,/P,, is calculated using 
the Same procedure as described earlier for Be eee Again, 
the value of Ee is checked on each iteration to ensure 
that it falls between O and 1 (line 780). If it does not 
then the interblade loss is reduced in a method analagous 
to that covered for the stator loss. 

After generating a value for P./Ps4 the program cal- 
culates GC: (line 1208). When the value for es is generated 
it is compared to 2 and three results are possible. 

i G. is equal to g, 2-05 then converganne has occurred and 
the program will continue. If @- is less than 6, (line 1250) 
then the value for interblade loss coefficient is increased 
by the current value of interblade increment and the program 
goes back to the start of the interblade calculations and 
begins anew. If eC. is greater than Cs (line 1240) the 

value of interblade loss is decreased in the same way 

aS was described for the stator loss. 

Two results are possible in the interblade region: 
Either the program finds the value of the interblade loss 
that gives flow angle convergence (@; = Gy)» or the inter- 


6) to 


blade loss increment becomes sufficiently small (10 
trigger termination. 

Note that the interblade loss increment is allowed to 
become very small (107°) (line 850). This is because the 


above mentioned angle convergence requirement is very 
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sensitive. Note, -also, that when the interblade loss 
increment falls to 107° and the iteration process in the 
interblade region is terminated, the stator loss coefficient 
is incremented (line 850). Stator loss is incremented 
because a higher value of stator loss may result in con- 
ditions that allow convergence in the interblade region. 

If not, the stator loss calculations eventually reach the 
termination criterion mentioned earlier. 

Once angle convergence has been attained in the inter- 
blade region then the value. of the rotor flow coefficient 
can be calculated (line 1470). If the program gets to this 
point for any given set of conditions then a solution 
exists and the only determination to be made is whether 
the resulting rotor flow coefficient is greater than the 
"choking" value (line 1530). If "ves" then it is dis- 
allowed, since a flow coefficient greater than the “choking” 
value cannot occur physically. If the calculated coef- 
ficient is less than, or equal to the "choking” value then 
the results are printed, the stator loss is incremented, and 
the calculations start again at the stator exit plane. 

The choking value of the flow coefficient for the 
rotor (2*, ) is determined in exactly the same way as the 
choking value for the stator (@*), with the assumed value 
of rotor loss used as the loss coefficient. (At present 
the value for ps is calculated separately and inserted in 
line 121. It is suggested that a "GOSUB 2140" be inserted 


at this point. Then the input value for rotor loss 
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coefficient can be changed easily and the "choking" value 
for the rotor will be calculated in the program.) If 

the resulting flow coefficient is greater than "choking", 
the program increments stator loss coefficient and returns 
to the stator plane calculations to begin anew. 

Therefore, in all cases the program will seek the 
highest value of stator loss coefficient for which solutions 
exist. At some point, the rotor flow coefficient normally 
attains a value less than the "choking" value and all 
solutions for higher stator loss coefficients will be 
aceeptable solutions. Or, if the "choking" condition is 
never reached, the program will attain the highest possible 
Stator loss in ae ene to achieve solutions less than 
"choking". 

When the highest value of the stator loss is reached, 
there are two possible avenues in the program. If the 
program was generating solutions less than "choking" at 
the time of termination in the stator calculations, then 
Kos will be incremented, stator loss will be set back to 
zero, and the process begun anew (line 2280). If the 
solutions being generated at the time of termination were 


greater than "choking", then the upper value of Ky has 


Ht 
been reached and any increase in Ks Will yield no further 
solutions. In this case Ka is incremented, Ks is set 

to the pre-determined lower limit, and stator loss is 


reset to zero (line 2330). 


The program continues in the above manner until the 
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pre-determined upper limit of k,_, is reached (line 2335), 


bi 
at which point the program is terminated. 

Since the 21-MX computer has only 6 digits of AcCCuUnaew, 
round-off error has been a problem. Where Newton's Method 
is used to solve for pressure ratio the computer lacks 
Sarticient accuracy romeenvenee to the desired accuracy 
under certain conditions. When this happens, the program 
will stay inside the iteration and never converge. This 
problem has been solved by putting a counter inside each 
iteration loop (lines 540 and 980). More than 60 iterations 
is treated aS a non-convergent condition, and the approp- 


riate loss coefficient is reduced in the manner covered 


previously. 


A-9 OPERATING PROCEDURE 

The procedures for operating the Hewlett-Packard 
21-MX computer will not be covered in this paper since 
the applicable manuals are available at the laboratory. 
The computer must be on, with the RTE-B operating system 
"READY". 

1. Load the paper tape program labeled "TTR 11". 

2. Edit into the program the’ minimum values of k 


bl 


(line 284) and k line (286). Edit the desired value of 


pa: 
AO, into line 280. If it is desired to get a print-out 

of every solution then either remove lines 1724 and 

1726 or place "REM" in front of them. If these two lines 


are left in the program then a solution will be printed 


out only if the calculated pressure ratio (Po/P..) is less 
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than, or greater than, the previously calculated lowest or 
highest pressure ratio respectively, for the value of Kd 
then being used. 

3. Type “RUN", “RETURN” on the keyboard and the 
program will begin execution by asking for inputs: 

"INPUT STATOR LOSS COEFFICIENT” The normal input is 
zero. However, if some specific case iS required then any 
desired loss coefficient can be input. 

"INPUT PTO” Input the upstream tetal pressure. 

"INPUT P2" Input the hood pressure. 

"INFUT RPM" Input the RPM desired. 

"INPUT TTO"” Input the total temperature into the 

Stage. 
This is all that is required to operate the program. 


Depending on the range of Ky requested the program will 


1 
take from 30 minutes to 36 hours to run on the 21-MX. 

It is important to realize that the above procedure 
will produce solutions that give values of the rotor 
flow coefficient less than "choking". If the turbine 
performance at the "choking" céndition is desired then 
after the computer run is complete a further step must be 
taken. It is necessary to take the values of Ky and Ks 
for which a solution is desired and force the choking 
condition. This is done as follows: 


1. Scan the output results at the desired k and 


bl 
Kyi and locate the point that has a value of rotor flow 


coefficient closest to the choking value and a non-zero 
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stator loss coefficient. 

2. Re-start the program, this time putting in a value 
Gretne Stator loss coefficient Slightly less than the value 
printed on the output. 

As stator loss is manually decreased, the rotor flow 
coefficient will increase towards the choking value. When 
the rotor flow coefficient is within the limits specified 
in line 2370 the words "CONVERGENT CONDITIONS" will appear 
on the atioitine output preceding the printed results. 

If the stator loss coefficient 1s decreased too much, 
the rotor flow coefficient will increase to a value above 
choking. In this case no ouput will be printed on the tele- 
type. However, the flow coefficient will appear on the 
video display. In this case increase the stator loss 
slightly. 

It has been found in this work that the extreme values 
of referred horsepower, stator loss, rotor loss, and 
efficiency occur at the lowest and highest values of KG 
for each value of Kt? With this in mind it is possible to 
find the range of predicted values at choking by forcing 
to choking only the lowest value of Ks (with non-zero 


Tor each. 


stator loss) and the highest value of Ky s bl 
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TABLE A-I 


DEFINITION OF VARIABLES 


LA*) Area at the Stator Throat 

(A, ) Area at the Stator Exit Plane 
Decision Variable 

Decision Variable 


Decision Variable 


Clee 


d (7) 


dF 
StF 


ref 


(AT; 5) 
Decision Variable 

(mg) Total-to-Static Efficiency 

(Q*) Flow Coefficient at Stator Throat 

(2, ) Flow Coefficient at Stator Exit Plane 
(2.) Flow Coefficient Through Rotor Passage 
(2. ) Flow Coefficient Through Interblade Area 
Input Value of 2. Based on Assumed Rotor Loss 
Decision Variable 

Decision V-riable 

(Gamma - 1)/Gamma 

Horsepower 

Referred Horsepower 

(§) Rho 

Reynolds Number 


98 





K1 
K4 
K5 
K6 
M1 
M2 


NO 
On 
03 
O4 


PO 
Pi 
ie 
P3 
P4 
P5 
P6 
P? 
P8 
PY 


R7 
R8 
RQ 


iz) eae 

Ks Blockage Factor at Stator Throat 

Koi Interblade Blockage Factor 

KOR Blockage Factor in the Rotor 

Ky Blockage Factor at the Stator Exit Plane 
Rotor Torque 

Referred Rotor Torque 

RPM 

Referred RPM 

(X,) Flow Angle at Stator Exit Plane 

(A;) Flow Angle at Rotor Entrance Plane 
(Ak; ) Increment for k, 
Pis/P to 

E 


1 


to 

Ps/P 46 

Po/P 45 

Py/P ey 

Pt1/P to 

Pi/P a4 

Pi/P eo 

Pai /P to 

Decision Variable 
26 
(feo) Total Density at Stagnation Conditions 
Relaxation Parameter 


Relaxation Parameter 


Predicted P./P,., 


ae 





S7-S9 Used in GOSUB 7000 to Find P,/Pp; 


TO 
Ea 
a2 
T3 
TY 
T5 
T6 
U1 
U2 
U3 
VO 
WO 
wi 
W9 
X0 
X1 
X2 
X3 
X4 
X5 
X6 
X8 
2 
Y6 
ae 
Y8-Y9 


TL0 

Teo/T to 
T/T 46 
T/T Lo 
Ta3/T to 
T/T 46 
T51/T to 


CG 


, Dimensionless Velocity (See Appendix A) 


> Dimensionless Velocity (See Appendix A) 


CQ 


U, Dimensionless Velocity (See Appendix A) 
on Total Velocity 


@ 


w Flow Rate 
w* Referred Flow Rate 


Decision Variable 


Xa Non-Dimensional Relative Velocity 
X, Non-Dimensional Velocity 
X, Non-Dimensional Velocity 
X. Non-Dimensional Velocity 
Xs Non-Dimensional Velocity 
X Non-Dimensional Relative Velocity 


Wi 
Xo Non-Dimensional Relative Velocity 
Decision Variable 

Non-Dimensional Velocity 

Increment for Ky 


Increment for Kos 


Counters 
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ZO 
Z3 
27 
Z8 
Z9 


Zo Stator Loss Coefficient 

Zs Interblade Loss Coefficient 
Z,_>5 Predicted Rotor Loss Coefficient 
Increment for Stator Loss Coefficient 


Increment for Interblade Loss Coefficient 
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PRESSURE DISTRIBUTION THROUGH THE STATOR 


FIGURE A-1i 
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FIGURE A-4 TYPICAL SOLUTION PROCESS 
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FIG. A-5 PROGRAM SCHEMATIC 
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APPENDIX B 


TURBINE TEST RIG (TTR) DATA REDUCTION AND PROCESSING 


B-1 INTRODUCTION 

This appendix describes only those changes that have 
been made to the data reduction procedure given in Ref. 2. 
The program numbers and the functions of those programs 
given in Ref. 2 have not changed. The channel and port 
assignments for data collection and storage have not 
changed. Variables also have not changed, although 
additional variables were pened in the course of mod- 
ifying the data reduction process. Those variables that 
now exist in addition to those given in Ref. 2 are given 
in Table B-I. 

Modifications made during the course of this work were: 

1. Raw data can now be read directly from mass memory. 
Previously it was necessary to read the paper tape for a 
given point in order to reduce it. 

2. Provisions have been written into the programs 
for smoothing the reduced data. 

3. All of the parameters used in calculating the 
theoretical loss coefficients that were previously 
entered from charts are now in polynomial form in the 
program. 

4, Storage of raw data is now a separate program. 
Previously it was necessary to run the reduction sequence 


to store data. 


iz 





B-2 DESCRIPTION OF PROGRAMS 

The data reduction is divided into ten separate 
programs with an additional program to store the raw data. 
Fig. B-1 shows the contents of each program and also shows 
the sequence in which the programs are chained. Following 
is a description of the reduction sequence referred to 
Fig. B-1 which will be followed by a description of how to 
run the reduction program. 

TTR in Fig. B-1 is used only to store the raw data. 
TTR does not chain to any other program. 

The data reduction process begins in TTR1. In order 
to smooth the reduced data as discussed in Section III it 
is necessary to have n values for (Pa-Prup)/(Peip Phun?’ 
Therefore, the first step in the data reduction process is 
to go from TTR1 to the point in TTR2 where the values of 
Ee, FP 


and P have been calculated. On completion, 


a’ hub’ 
the program returns to TTR1 and repeats the process n 


tap 


times until the n values of (P,-P ari Ee Bae) have 
been calculated. The program then proceeds to TTR9 where 
a polynomial curve fit for the n values of 


(Ps -Prup Pian? hu 
returns the polynomial coefficients to TTR1B where the 


») is generated. Then the program 


reduction process begins at the first point in the run. 
However, when the point in TTR2 is reached where Pe is 
calculated, it is calculated using the polynomial as 
described in Section III. 


The rest of the reduction process in Fig. B-1 is the 


Iv 





Same aS described in Ref. 2 with the exception that no 
further keyboard inputs are needed until the raw and 


reduced data is tabulated. 


B-3 RAW DATA STORAGE 

The procedure for storing raw data on the mass memory 
is outlined in Ref. 2, Appendix B. Sec. B-3, lines 10719. 
Note that line 10 should read, "GET ‘TTR‘'". Although it is 
now a Separate program the procedures for storing data 


have not changed. 


B-4 PROCEDURE FOR DATA REDUCTION PROGRAM 
1. Key in GET "TTR1B* 
2. Press “RUN EXECUTE”. 
3. The following check list will be printed on the 
HP 9830 printer: 
"PRIOR TO RUNNING THIS SEQUENCE ENSURE FOLLOWING: 
"THAT TIR2 LINE, 590g14S PRORER FACTOR IN IT; 
1.02 FOR HOODED 
1.01 FOR UNHOODED 
"IF IT IS NECESSARY TO INPUT ALPHA 1 THEN CHANGE TTR2 
"LINE 960 TO ‘INPUT A3° 
"IF BLOCKAGE FACTOR OTHER THAN 1, CHANGE TTR2-1060 
"TQ ‘INPUT X7°, AND PUT SEMICOLON AFTER TTR2-1050 
"ENSURE THAT TTR1, TTRiB HAVE THE PROPER FILENAME IN 441 
"IF IT ISN'T DESIRED TO STORE REDUCED DATA CHANGE 
TTR6-580 TO 'G1=0'° 
“ENSURE TTR6-610, 630 HAVE PROPER FILENAME FOR REDUCED 


DATA 
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"ENSURE TTR7-220 HAS PROPER FILENAME FOR RAW DATA 

"ENSURE TTR8-220 HAS PROPER FILENAME FOR REDUCED DATA 

"IF IT IS DESIRED TO RUN A SINCE PCrNT. = wit noun, 

SMOOTHING, INSERT TTR2-1127 ‘GOTO 1200°. OTHERWISE 
OMIT Pik 2127 

The above checklist will prepare all ten programs. The 
last item is a provision for elimination of smoothing if 
it is desired to run a single point. If smoothing is not 
desired then each point must be run separately with 
TTR2-1127 inserted. 

4, Press CONTINUE EXECUTE 

GS, The display will read LOWEST, HIGHEST RECORD # 

THIS RUN. Input the lowest record number and the 
highest record number from which it is desired to read 
raw data. Note, the record number on which the reduced 
data will be stored will be the same number as the raw 
data record from which the raw data was read. Make sure 
the filename designated for the storage of the reduced 
data can accept data on those record numbers without 
writing over previous information. 

6. Next the display will read PRINT OUT RESULTS? 
YES=1, NO=0O. If ailis selected here then at the conclusion 
of data reduction, for all the points in the run, the 
program will begin to print all results. If 0 is selected 
then all reduced data will be stored on the mass memory. 

7- It will take about 15 minutes for the program to 


read the raw data and make the calculations necessary for 


iD 





smoothing the data. When this is oompleted the program will 
start reducing each point and printing the results. 

8. At the completion of data reduction for all 
points the option to tabulate the data is available. 
ENTER RECORD #'S: LOWEST, HIGHEST will appear on the display. 
If it is desired to tabulate the raw data then enter the same 
record numbers that were entered to initiate the program. 
The raw data will then be tabulated. 

9. Next the display will read ENTER LOWEST, HIGHEST 
RECORD NUMBER. If it is desired to tabulate the reduced 
data then enter the appropriate record numbers. The 


reduced data will be tabulated. 
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Az 


A5 


A6 


A? 


F8 


TY 


TABLE B-I 
VARIABLES ADDED TO DATA REDUCTION PROGRAM 


Lowest Record Number 

Counter 

Peso ne ne 

Highest Record Number 

(Array) Polynomial Coefficients From TTR9 
Decision Variable 

(Array) Values of Isentropic Head Coefficient 
Polynomial Approximation of Reynolds Number 


Decision Variable 


Ale, 





isd 


1.) RAW DATA STORAGE 


2.) RAW DATA PRINTED 









TTR 9 
1.) GENERATE POLYNOMIAL 


CURV Sw =f OR 1 


Pin-Py 







VALUES OF 






TTR 1B 
1.) READ RAW DATA 





1.) READ RAW DATA 






2.) CHANNEL ASSIGNMENT 2.) CHANNEL ASSIGNMENT 


3.) TEMPERATURE SUBROUTINE 3.) TEMPERATURE SUBROUTINE 





TR 2 
1.) CALCULATE MASS 


FLOW RATE 


2.) EVALUATION OF 
CONTROL VOLUME A 


ERG 
1.) EVALUATION OF CONTROL 


VOLUME B AND TEMPERATURES 





FIGURE B-1 DATA REDUCTION SCHEMATIC 


118 





TTR 4 
1.) CALCULATE STATOR 


EA een Galea 
2.) CALCULATE THEORETICA 
LOSSES 





es 
1.)PRINT REDUCED DATA 


TTR 6 
1.) STORE REDUCED DATA 


1.) TABULATE RAW DATA 


TTR eG 


1.) TABULATE REDUCED DAT: 
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